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LONG-TERM  GOALS 

This  is  the  final  year  of  the  four  year  grant.  My  main  objectives  in  this  final  year  have  been  to  carry 
out  the  following  studies:  to  publish  the  findings  of  the  now  completed  work  on  finding  reliable 
estimates  of  the  theoretical  limits  to  tropical  cyclone  forecast  errors,  for  both  track  and  intensity 
predictions,  but  with  emphasis  on  distinguishing  between  overall  track  errors,  and  position  and  timing 
errors  at  landfall;  to  complete  the  work  that  has  addressed  the  problem  of  identifying,  quantifying  and 
reducing  these  sources  of  errors;  to  continue  to  improve  the  initial  state  specification  for  TCs,  by  the 
development  and  application  of  4D  data  assimilation  procedures,  particularly  with  respect  to  the 
ingestion  of  newly  emerging  data  sources,  especially  satellite,  radar  and  scatterometer  data;  to  develop 
and  implement  major  model  changes  and  additions,  motivated  by  the  very  active  2004  season;in  the 
Atlantic  Ocean,  the  Caribbean,  and  the  Gulf  of  Mexico;  to  complete  the  work  on  comparing  existing 
and  new  methods  for  the  evaluation  of  the  skill  of  TC  track  prediction  schemes,  part  of  which  is  to 
develop  new  probabilistic  statistical  and  statistical-dynamical  procedures  for  TC  forecasting  beyond 
24  hours,  and  also  to  include  second  order  and  non-linear  effects;  and,  finally,  to  complete  work  on  a 
specific  impact  of  landfalling  TCs,  namely,  the  generation  of  tornadoes,  especially  tornado  outbreaks. 
My  goals  all  have  implications  for  transitions. 

OBJECTIVES 

As  mentioned  in  previous  section,  there  are  six  main  scientific  objectives  of  this  proposal.  My  first 
objective  is  to  publish  our  findings  on  estimating  the  intrinsic  limits  of  predictability  of  TC  cyclone 
behavior,  for  both  mean  absolute  forecast  position  errors  and  intensity  errors.  As  mentioned  in 
previous  FY  reports,  intrinsic  limits  of  predictability  for  TCs  exist  because  the  equations  governing  the 
behavior  of  atmospheric  systems,  including  TCs,  are  deterministically  chaotic.  Errors  in  initial 
conditions,  model  formulation,  and  boundary  conditions  generate  errors  that  grow  until  the  predictions 
have  no  skill.  The  first  objective  leads  to  the  second,  which  is  to  carry  out  research  directed  at 
identifying  and  understanding  the  sources  of  errors  in  both  the  initial  conditions  and  the  model 
formulation,  and  to  reducing  these  errors.  This  objective  remains  the  major  part  of  the  program  of 
numerical  analysis  and  prediction  of  this  proposal.  My  focus  largely  will  be  on  TCs  that  are 
approaching  or  making  landfall,  as  they  are  the  most  devastating  in  terms  of  loss  of  life  and  property. 
My  third  objective  is  to  continuing  the  development  of  my  4D-VAR  data  assimilation  scheme.  The 
sub-components  of  the  third  objective  continue  to  be  to  improve  upon  the  present  4D-VAR  scheme, 
especially  its  efficiency,  and  its  ability  to  ingest  the  many  emerging  data  sources  such  as  the  satellite 
and  radar  data.  My  fourth  objective  has  been  to  complete  an  evaluation  of  existing  schemes  for 
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assessing  NWP  model  performance,  and  to  develop  alternative  procedures.  Present  methods  of 
assessment  frequently  are  either  biased  towards  the  NWP  models  or  are  too  scale  dependent,  unduly 
rewarding  smooth  forecasts  and  penalizing  detailed  predictions.  In  particular,  new  assessment 
schemes  are  required  for  TC  model  forecasts,  as  they  are  not  well  handled  by  conventional  measures  of 
skill,  given  the  many  scales  and  the  range  of  variables  that  need  to  be  assessed.  As  ensemble 
techniques  dominate  forecasting  systems,  the  need  for  new  assessment  procedures  has  increased.  A 
fifth  objective  is  the  continued  development  of  a  new  statistical  technique  for  either  stand-alone 
forecasting  or  to  correct  TC  model  predictions.  Finally,  a  sixth  objective  was  added  in  FY04,  and  has 
proved  to  be  timely  given  the  events  of  the  2004  hurricane  season.  This  objective  is  a  careful  analysis 
of  tornado  outbreaks  associated  with  landfalling  TCs. 

APPROACH 

My  approach  for  the  six  goals  is  as  follows.  The  methodology  employed  for  the  first  goal  is  explained 
in  the  ONR  FY01  and  FY02  reports.  Two  distinct  techniques  have  produced  almost  identical  answers. 
One  approach  was  to  generate  an  ensemble  of  initial  model  states  using  the  archived  data  sets  from 
various  operational  global  NWP  centers  around  the  world.  The  second  technique  was  based  purely  on 
observational  data  in  which  the  spread  rate  of  initially  close  pieces  of  TC  trajectories.  These  methods 
have  now  been  applied  to  the  problem  of  estimating  the  predictability  characteristics  of  TC  intensity 
and  intensity  change.  The  task  has  been  massive  computationally,  but  it  is  vitally  important  to  know 
how  close  we  are  in  practice  to  the  limits  of  predictability.  The  second  goal  has  been  to  obtain  model 
simulations  with  much  more  realistic  TC  structure,  intensity,  intensity  change  and  motion  as  TCs 
approach  landfall.  The  second  goal  links  directly  with  my  third  goal,  as  the  approach  in  both  cases  is 
use  my  HIRES  data  assimilation  and  NWP  model  systems.  In  FY04  I  continued  to  use  4-D  variational 
assimilation  procedures  without  TC  bogussing,  using  high  spatial-temporal  frequency  satellite  derived 
data  of  various  types  and  as  many  other  sources  of  data  as  possible,  especially  radar  data.  Populated 
coastal  regions  are  threatened  each  year  by  landfalling  TCs.  Therefore,  it  is  of  primary  importance  to 
acquire  a  better  understanding  of  their  internal  dynamics  in  order  to  improve  forecasting  capabilities. 
Here,  emphasis  is  placed  on  the  importance  of  improving  our  knowledge  of  the  TC’s  microphysical 
and  electrical  structure.  Further  motivation  for  this  work  comes  from  observational  studies,  which 
stressed  the  importance  of  more  systematic  monitoring  of  any  change  in  TC  lightning  flash  activity. 

The  latter  could  be  a  very  useful  source  to  detennine  the  distribution  of  the  convective  precipitation, 
particularly  at  landfall,  where  it  could  be  used  to  provide  advance  warnings  for  potential  flooding  and 
severe  weather  originating  from  embedded  supercells.  The  approach  to  the  fourth  goal  has  been  to 
compare  existing  measures  of  assessment  of  NWP  model  skill  with  alternative  measures.  Initially,  the 
alternative  measures  were  standard,  but  have  now  been  extended  well  beyond  the  existing  measures. 
The  approach  to  the  fifth  goal  has  been  to  use  the  Atlantic  basin  TC  data  base  to  generate  conditional 
probabilities  for  predicting  the  future  tracks  and  intensities  of  TCs  in  that  basin.  An  example  is  given 
in  the  Results  section  below.  Finally,  the  sixth  goal  is  to  define  and  develop  a  tornado  outbreak 
climatology  from  1953  to  2002  over  the  US  for  tornadoes  associated  with  land-falling  hurricanes. 
Tornadoes  generated  by  TCs  tend  to  be  weaker  and  land-falling  hurricanes  can  persist  for  several  days. 
Finally,  case  studies  that  meet  the  outbreak  definition  compared  against  several  non-hurricane  related 
tornado  outbreaks  in  the  U.S.  TC  Beulah  (1967)  holds  the  current  record  with  1 17  tornadoes  reported. 
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WORK  COMPLETED 


My  work  on  the  first  goal  is  now  completed  and  the  work  has  been  published.  Recent  work  has 
focussed  on  reducing  mean  track  errors  for  TCs  that  are  regarded  as  difficult  to  predict.  Applying  the 
same  procedure  to  TC  intensity  predictability  has  also  been  completed,  with  the  emphasis  again  being 
on  the  more  difficult  storms.  This  work  has  produced  some  intriguing  and  unexpected  results,  most 
notably  the  large  variations  from  basin  to  basin.  Second,  development  of  the  HIRES  data  assimilation 
and  prediction  system  has  continued,  with  the  incorporation  of  new  data  from  radar  and  satellite 
sources.  This  work  also  has  been  published.  My  third  goal  of  understanding  the  factors  contributing  to 
improving  the  forecasts  of  TC  intensity  and  structure  through  improvements  to  the  NWP  model  have 
reached  a  new  phase  with  the  successful  introduction  of  an  improved  cloud  microphysics  scheme.  The 
inclusion  of  high  quality  cloud  microphysics  in  the  NWP  models  appears  to  be  vital  for  understanding 
and  predicting  TC  behaviour.  This  is  true  especially  as  the  TCs  approach  landfall.  Some  significant 
findings  have  been  made  and  will  be  presented  at  the  2005  AMS  annual  meeting.  Fourth,  existing 
methods  for  assessing  NWP  model  performance  have  been  shown  to  be  far  too  simplistic  and  also 
must  be  extended  to  cover  ensemble  predictions,  which  are  now  becoming  common-place  in  usage. 
New  techniques  have  been  tested  and  results  also  will  be  presented  at  the  2005  AMS  meeting. 


Hurricane  Ivan 
September  9, 2004 
11  AM  EDT  Thursday 

NWS  TPC/National  Hurricane  Center 
Advisory  29 

Current  Center  Location  14.5  N  71.4  W 
Max  Sustained  Wind  160  mph 
Current  Movement  WNW  at  15  mph 
(S)  Current  Center  Location 
•  Forecast  Center  Positions 
H  Sustained  wind  =■  73  mph 
L2~~\  Potential  Day  1-3  Track  Area 
H  Hurricane  Warning 
Hurricane  Watch 
Hi  Tropical  Storm  Warning 
Tropical  Storm  Watch 


Tropical  Cyclone  Motion  Climatology 
from  0  to  24  hrs 

gray  circles  represent  frequency  per  sector 


Average  Cyclone  Speeds:  n  m 

0  to  5  kts 
6  to  10  kts 

11  to  15  kts 

>20  kts 

Latitude:  16 . 5N  Longitude:  75. 1W  Radius:  500. km 
Climatology  filtered  by  location.  n=  134 


Figure  1:  The  official  track  forecast  from  NHC  (left),  compared  with  the  probabilistic  prediction 
(right).  The  probabilistic  prediction  suggested  a  more  westerly  movement  of  Hurricane  Ivan, 

which  did  in  fact  occur. 

Fifth,  the  probabilistic-statistical  prediction  scheme  has  been  tested  on  every  stonn  throughout  the 
2004  hurricane  season  in  the  Atlantic  basin.  It  has  been  particularly  useful  in  forecasting  short  term 
track  directions  and  speeds  of  movement.  In  a  number  of  cases,  the  scheme  gave  quite  different 
predictions  from  the  NWP  models  and  the  official  NHC  forecast.  On  every  occasion,  the  probabilistic 
scheme  was  more  accurate  over  the  first  24  hours.  A  journal  paper  is  in  preparation  on  the  performance 
of  the  scheme  over  the  2004  season,  and  also  will  be  presented  at  the  2005  AMS  annual  meeting.  An 
example  is  shown  above  (in  Figure  1),  for  Hurricane  Ivan,  when  located  at  14. 5N,  71.4W.  The  official 
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forecast  moved  Ivan  northwest.  The  statistical  scheme,  based  on  134  storms  moved  Ivan  due  west  on  a 
path  that  had  it  entering  the  Gulf  of  Mexico,  well  to  the  west  of  Florida.  This  is  indeed  what  occurred. 

RESULTS 

Results  are  now  complete  for  the  first  goal,  which  was  to  determine  how  close  current  NWP  models 
are  to  our  best  estimates  of  the  limits  of  predictability  for  mean  absolute  TC  track  errors.  The  major 
finding  of  FY01,  FY02  and  FY03  was  the  potential  for  TC  track  errors  to  at  least  be  halved.  A  second 
set  of  results  concerns  the  performance  of  the  NWP  systems  presently  at  the  core  of  my  ONR  research 
program.  This  model  is  the  High  Resolution  (HIRES)  data  assimilation  and  prediction  system. 

Largely  using  ONR  support,  this  system  has  been  developed  over  six  years,  for  use  in  TC  applications 
by  the  PI,  various  graduate  students  and  part-time  time  research  associates.  It  has  continued  to 
undergo  development  in  FY04,  especially  in  the  area  of  data  assimilation  and  representation  of  cloud 
microphysics. 


Fig.  2.  Time-height  plot  of pressure  at  z  =  100  m  AGL  for  the  first  27.5  hours  of  simulation.  Blue 

(red)  line  shown  for  the  land  (no-land)  case. 
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At  26h  40  min  hours  (1600  min)  of  simulation,  the  vortex  shows  a  central  pressure  of  946  mb  (Fig.  2, 
redline)  and  a  distinct  eyewall  with  several  connecting  rainbands  (Fig.  3a).  At  that  time,  tangential 
winds  are  exceeding  50  m  s'1,  with  the  highest  gusts  just  within  the  boundary  layer  above  ground, 
corresponding  to  a  low-end  Category  3  storm.  Consistent  with  many  in  situ  observations  the  TC 
eyewall  exhibits  an  outward  tilt  in  reflectivity  and  in  the  tangential  wind  profile.  As  the  storm  weakens 
over  land  tangential  wind  profile  becomes  more  vertically  stacked.  The  weakening  of  the  TC  is  clearly 
evident  for  instance  in  the  central  pressure  at  100m  above  the  ground  (Fig.  2,  blue  line)  and  in  the 
tangential  wind  profiles  (Figs.  3c  and  3d),  where  maximum  meridional  (tangential)  wind  component 
are  exceeding  50  m  s'1  over  the  ocean,  while  remaining  below  40  m  s'1  once  over  land  at  the  same 
time.  The  horizontal  reflectivity  field  of  the  TC  over  land  also  show  many  interesting  features 
consistent  with  in  situ  observations,  namely,  weaker  reflectivity  in  and  around  the  eyewall  and  the 
absence  of  a  closed  circular  eyewall  and  eye  near  and  around  the  storm’s  center 
(compare  Figs.  3a  and  3b).  This  is  a  consequence  of  weaker  updraft  mass  flux  in  the  eyewall  and  hence 
weaker  divergence  aloft,  in  turn  reducing  the  subsidence  and  the  subsequent  warming  responsible  for 
the  formation  of  the  eye.  On  the  other  hand,  when  the  storm  evolves  over  open  waters,  more 
continuous  convection  in  the  eyewall  and  a  well-defined  closed  circular  eye  are  noticeable  near  the 
storm’s  center.  Between  these  two  convective  regions,  a  stratiform  region  with  weak  reflectivity  values 
is  also  evident.  Secondary  weak  bands  underneath  the  outer  core  stratiform  region  of  the  stonn’s 
eyewall  are  also  found  on  the  storm’s  western  and  eastern  flank. 
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Fig.  3.  (a)  Horizontal  cross  section  of  radar  reflectivity  (dBZ)  shown  by  5  dBZ  contour  increments 
from  5  dBZ  to  75  dBZ  at  t  =  1600  min  for  the  no  land  case,  and  (b)  the  land  case,  at  z  =  1,017  km. 
Locations  of  CG  lightning  strikes  are  also  shown  by  a  cross  for  -CG  flashes  and  by  a  +  symbol  for 
+CG  flashes.  The  flash  locations  were  plotted  for  6000  s  time  interval  from  the  cross  section  time. 
Panels  (c)  and  (d)  show  vertical  cross  section  in  the  X-Z  direction  of  the  V  (or  Y)  component  of  the 
wind  field  aty  =  330  km,  with  corresponding  scale  shown  below. 
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IMPACT/APPLICATIONS 


In  FY04,  the  final  phase  of  the  estimation  of  how  close  we  are  to  the  limits  of  track  predictability  on 
the  0-4  day  time  scale  have  now  confirmed  earlier  work  that:  large  TC  track  forecast  errors  remain;  and 
that  these  track  errors  are  still  at  least  twice  that  which  can  theoretically  be  achieved.  The  patchy 
performance  of  NWP  models  over  the  2004  season  appears  to  have  confirmed  these  findings.  The  gap 
between  that  still  exists  between  practice  and  theoretical  limits  points  to  the  necessity  for  continued 
understanding  of  TC  processes  and  improvements  in  the  analyses  and  NWP  model  formulation.  A 
second  impact  in  FY 04  resulted  from  the  continued  use  of  careful  quality  controlled,  high  resolution, 
data  from  existing  and  new  sources;  developments  in  continuous  (4D)  assimilation;  and  ongoing  model 
improvements  and  increased  model  resolutions  of  5km  or  below.  TC  simulations  that  have  many  of 
the  observed  features  of  actual  TCs  will  soon  be  routine,  as  suggested  above  in  Figure  3.  Also,  the 
adjoint  sensitivity  techniques  developed  earlier  continue  to  improve  our  ability  to  predict  TC  motion, 
structure,  intensity,  and  intensity  change.  Fifth,  the  range  of  precursors  of  TC  activity  and  intensity 
identified  in  FY02  in  the  eastern  Indian  Ocean,  northwest  Australia  and  the  western  Indian  Ocean  have 
been  extended  to  the  western  North  Pacific  Ocean  basin.  These  precursors  were  again  shown  to  be 
present  at  the  inter-seasonal  level,  through  the  intra-seasonal  level,  to  time-scales  of  days  and  even 
during  the  life  cycles  of  TCs.  Other  applications  include  model  improvements  which  should  increase 
our  understanding  of  changes  in  TC  structure,  and  associated  severe  weather,  as  they  approach  and 
make  landfall;  and  a  return  to  statistical  schemes  that  provide  rapid  updating  of  probabilistic  forecasts 
of  TC  track  speed  and  direction. 

RELATED  PROJECTS 

I  am  a  co-PI  on  the  CBLAST  grant  entitled  “The  Impact  of  Air-Sea  Interaction  Research  on  Larger- 
Scale  Geophysical  Flows.”  with  Dr  Michael  Banner  of  The  University  of  New  South  Wales.  I  also 
continue  to  interact  with  Dr  Wang  Yuqing,  University  of  Hawaii,  on  the  impact  of  cloud  microphysics 
schemes. 

SUMMARY 

This  proposal  aims  at  increasing  our  understanding  of,  and  ability  to  predict,  tropical  cyclones  in  a 
number  of  areas.  Tropical  cyclones,  which  are  also  referred  to  as  hurricanes  and  typhoons,  are  the 
most  devastating  storms  on  earth.  As  such  it  is  vital  to  understand  and  to  predict  their  behavior.  To  do 
so  in  an  accurate  and  timely  manner  requires  research  on  their  motion,  their  structure,  their  intensity, 
especially  when  nearing  land.  To  achieve  these  aims,  a  program  of  data  collection  and  computer 
model  simulations  is  being  carried  out,  with  the  predictions  being  compared  with  observations  of 
selected  storms.  Deficiencies  in  the  initial  data,  the  model  formulation  and  the  model  predictions  are 
then  identified  and  research  is  carried  out  on  improving  these  aspects.  The  ultimate  goal,  expressed 
succinctly  is  accurate,  timely  and  reliable  model  predictions  of  TC  tracks  and  intensity,  especially  for 
storms  that  threaten  coastlines.  This  ONR  funded  research  program  will  help  make  the  School  of 
Meteorology,  University  of  Oklahoma,  an  international  focal  point  in  tropical  cyclone  research. 
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